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This study is devoted to the experimental analysis of the Strato-rotational Instability (SRI). 
This instability affects the classical cylindrical Couette flow when the fluid is stably stratified in 
the axial direction. In agreement with recent theoretical and numerical analyses, we describe for 
the first time in detail the destabilization of the stratified flow below the Rayleigh line (i.e. the 
stability threshold without stratification). We confirm that the unstable modes of the SRI are non 
axisymmetric, oscillatory, and take place as soon as the azimuthal linear velocity decreases along 
the radial direction. This new instability is relevant for accretion disks. 
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Taylor- Couette flow (or cylindrical Couette flow) is cer- 
tainly one of the most popular laboratory flows and its 
study has already led to an abundant scientific literature. 
Several articles review the theoretical and experimental 
features of this flow that consists of a simple shear be- 
tween two co-axial rotating cylinders [1]. Its linear sta- 
bility is given by the Rayleigh criterion which states that 
rotating shear flows are stable when their angular mo- 
mentum increases radially. The richness of the transi- 
tion diagram that reports the different flow patterns and 
their subsequent transitions to chaos and turbulence has 
been extensively explored by Andereck et al. Amaz- 
ingly, the literature about Taylor- Couette flows in the 
presence of an axial stable stratification is much more 
limited even though it finds direct applications in geo- 
physics (atmospheric or oceanic flows) or astrophysics 
(accretion disks for instance). The first experiment on 
stratified Taylor- Couette flow was performed by Thorpe 
[3] in 1966. His results were further extended by the ex- 
periments of Withjack and Chen [4], Boubnov et al. [5|, 
Caton et al. [6] and by the numerical simulations of Hua 
et al. [7']. Most of these studies were performed with a 
stationary outer cylinder and all concluded that stable 
vertical stratification stabilizes the flow and reduces the 
axial wavelength of the unstable modes which bifurcate 
through Hopf bifurcations. 
In astrophysics, Keplerian flows whose angular velocity 



is Vt{r) 



-3/2 



, are always stable to infinitely small per- 



turbations in respect to the Rayleigh criterion [8]. As a 
consequence of this stability, several scenarios have been 
recently proposed to justify the existence of turbulence in 
accretion disk flows. Ingredients other than pure shear 
then need to be taken into account in order to trigger 
instabilities and transition to turbulence in these flows. 
One of the most popular scenarios relies on the Magneto- 
Rotational Instability [9] but other instabilities such as 
the elliptical instability [10] or non linear shear instabil- 
ities [8j have also been suggested. 

In 2001, Molemaker et al. [ll| and Yavneh et al. [H 



predicted that cylindrical Couette flows in a stratified 
fluid may become unstable even if the Rayleigh criterion 
for stability was verified, i.e. in the corresponding sta- 
ble regime of pure fluid flow. Moreover, and contrary to 
the classical Taylor vortices of the centrifugal instability, 
the most unstable modes should be non- axisymmetric. 
This theoretical analysis has then been continued in an 
astrophysical context by Shalybkov and Riidiger [13] and 
extended to the stability of accretion disk Keplerian flows 
by Dubrulle et al. [14] who gave to this new instability 
its now accepted appellation Strato-Rotational Instabil- 
ity (SRI). To the best of our know ledg e, and despite ex- 
plicit calls recently published [ll|, Il3l | no experimental 
validation of these theoretical predictions has yet been 
provided, apart from a short comment made by With- 
jack and Chen [3] in 1974 and a stability curve that was 
determined by Boubnov and Hopfinger in 1997 [lH. We 
present here the first quantitative experimental evidences 
of the Strato-rotational Instability and its expected heli- 
coidal modes. 

Our study is based on Kalliroscope visualizations per- 
formed in a classical Couette device with a ratio of inner 
to outer cylinder diameters equal to t] = 0.80. The shear 
created by the differential rotation of these cylinders is 
measured by the ratio fi of the outer cylinder rotation 
speed to the inner cylinder rotation speed: /i = 
Using /i and 7^, the two traditional control parameters 
for the cylindrical Couette flow, the Rayleigh criterion 
reads /i = which is represented by the Rayleigh line in 
the (Qo^^i) plane. A salt stratification is realized along 
the axis of rotation of the Couette apparatus using the 
classical "double-bucket" filling method [16]. The cor- 
responding Froude number of the fluid is chosen equal 
to 0.5 which is the estimated value for accretion disks 
[1^. Contrary to the extensive study of Caton et al. Q 
where the outer cylinder is at rest, the particularity of 
our study comes from the possibility to rotate separately 
both cylinders and therefore to study the development 
of non axisymmetric helicoidal modes in usually stable 
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FIG. 1: Kalliroscope visualization of the stable Couette flow 
in pure water (a) and of the SRI in the vertically stratified 
salt water with Fr = 0.5 (b) at the same Reynolds number 
= 1155 (flo = 1.07rad/s, Q^ = 1.5rad/s). The braid pat- 
tern of the SRI is formed by the superimposition of two iden- 
tical helicoidal waves (azimuthal wavenumber m = 4) propa- 
gating in opposite vertical directions. The laser sheet shows 
the meridional section of the travelling pattern (on the right 
part of the device) . 



regions (see fig. 1). As we will see, our experimental de- 
termination of the Strato-rotational Instability threshold 
is in excellent agreement with the value /i = ?7 proposed 
in [13] and corresponding to a constant azimuthal linear 
velocity along the radial direction. 

Our Couette device consists of two co-axial cylinders 
whose length is 168mm. The inner radius of the outer 
cylinder is Rq = 69mm and the radius of the inner cylin- 
der is Ri = 55mm, conferring a 14mm gap between the 
two concentric cylinders. The ratio of the inner to outer 
radius of the fluid cavity is equal to 77 = 0.8 explic- 
itly chosen to allow some comparisons with the results 



of Shalybkov and Riidiger [T^. The cylinders are po- 
sitioned vertically. The inner cylinder is made of black 
polished acetal plastic and the outer one is of transpar- 
ent glass so visualizations can be performed with the help 
of Kalliroscope flakes. To complement these bulk visu- 
alizations that use a classical light bulb for illumination, 
a laser sheet can also be installed to get a description 
of the hydrodynamical structure in a vertical meridional 
plane (see figure 1). The cylinders are driven by a d.c. 
servo-motor. Rotational speeds are measured by an op- 
tical encoder with an accuracy better than 1%. The top 
and bottom lids of the device are fixed to the external 
cylinder and rotate with it. The experimental protocol is 
the following. The classical double-bucket technique [l^ 
is used to obtain a vertical stratification of salt water. 
The vertical density gradient is measured before each run 
with the help of a density meter or a conductivity meter 
using a calibrated plunging probe. A linear gradient is 
quite easily achieved and its value leads to the determi- 
nation of the Brunt- Vaisala pulsation: N = (— ^ ff )^^^ 
with an accuracy better than 5 %. When changing the 
salt concentration in one of the filling bucket, N can be 
chosen between and 3.1rad/s. The Froude number Fr 
is then calculated by the ratio of the inner cylinder rota- 
tion speed Qi to N. We have fixed a priori its value to 
Fr = 0.5 in order to compare our experimental results 
with the numerical predictions of [13i]. Once the density 
gradient is established, the inner cylinder rotation speed 
Qi is thus fixed for each run, which also fixes the value 
of the Reynolds number: Re = (Ro-Ri) ^ where u is 

the viscosity of water. For each run, the initial condition 
is solid body rotation obtained when both cylinders ro- 
tate at the same speed, then we systematically change the 
angular velocity ratio /i = ^ by slowly decreasing Qq- 
When the threshold for the SRI is reached a single travel- 
ling helicoidal wave appears in the flow. After some min- 
utes, a vertical counter propagative wave, having identi- 
cal azimuthal, temporal and vertical (in absolute value) 
frequencies is superimposed on the first one, and creates 
a braid pattern as illustrated in figure 1-b). We believe 
that the initial dissymmetry between both waves comes 
from initial conditions and that the reflection on the lids 
of the first wave restore the symmetry between -\-k and 
— k waves expected from the theory. Figure 1-a) shows 
the pure water cylindrical Couette flow in the exact same 
conditions and demonstrates the stability of the flow in 
the absence of stratification. Let us note that whereas it 
is known that stratification damps vertical motions and 
instabilities in some cases 0, [gI, 0, [Hj , it is quite clear 
here that it can also destabilize shear flows. 
Ten series of experiments corresponding to ten different 
Reynolds numbers ranging from 339 to 1210 have been 
performed. Figure 2 sumarizes the experimental stabil- 
ity diagram of the flow in the parameter plane (/i,i?e)- 
It shows our experimental data points (x for instabil- 
ity, o for stable flow) superimposed on the stability dia- 
gram calculated in [13]. As can be observed an excellent 
agreement is obtained. In particular, we confirm that the 
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FIG. 2: Stability diagram of the Strato- Rotational Instabil- 
ity in the (Re, /J.) parameter plane for a Froude number of 0.5 
obtained for N ranging from 0.85rad/s to S.lrad/s. o and 
X stand respectively for experimental observations of stable 
and unstable flows. Squares stand for stationary axisymetric 
Taylor vortices. These experimental data points are superim- 
posed on the solid lines of the theoretical stability diagram 
of Shalybkov and Riidiger [13] for azimuthal wavenumbers 
m = (e.g. the classical Taylor vortices) and m = 1,2,3,5 
for ?7 = 0.78 (some helicoidal modes of the SRI). Vertical lines 
correspond to the theoretical inviscid limits with respectively 
the Rayleigh line fi — rf relevant for the unstratified case 
(solid line) and the threshold for the stratified case /a = r] as 
suggested in (dotted line). 



inviscid threshold, proposed by Shalybkov and Riidiger 
pisj, is given by /i = 77. Further decreasing /i, patterns be- 
come increasingly complex as additional helicoidal modes 
(m 7^ 0) are progressively added. We suggest that the he- 
licoidal modes calculate d by Molemaker et al. [H] and 
Shalybkov and Riidiger [1^ and observed in our exper- 
iments should be by continuity the "vortex modes" de- 
scribed by Caton et al. [6] when the outer cylinder does 
not rotate. We never found the "standing wave" regime 
of Caton et al. |[6i] (i.e. their pulsating axisymmetric 
modes) within the explored range, presumably because 
its unstable region is too narrow and our experimental 
device is not precise enough to excite this very marginal 
state (see the comparison between the results of [s] and 
[6]). Nevertheless, in experiments at /i = and large 
Reynolds number, we observed the reappearance of sta- 
tionary axisymmetric Taylor vortices in agreement with 
the regime diagram of Boubnov et al. [5] . This stationary 
axisymmetric mode was never excited when decreasing /i 
for /i > 0, but once excited at /i = 0, they could then 
be maintained up to /i = 0.2. These interesting non- 
linear effects are beyond the scope of the present paper 
devoted to the SRI and will be further studied. Figure 3- 
a) shows the space-time diagram of the braid pattern ob- 
served at SRI threshold and shown in figure 1-b). These 
diagrams are built by assembling vertical lines from suc- 
cessive video images, and allow to extract in a systematic 
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FIG. 3: (a) Space- time diagram at threshold of the growing 
helicoidal mode presented in figure 1-b), with Fr = 0.5 and 
Re = 1155. (b) 2D-Fourier Power spectra of the space-time 
diagram: contours are equally space from to the maximum 
energy ^^"^ with a spacing of ^'^'^^/20. Non Unear interac- 
tion between both helicoidal waves is visible as peaks around 
^:^0 and 1.6s-\ 

way the temporal frequency / and the axial frequency q 
of the hydrodynamic structures using a 2D-Fourier trans- 
form, as shown in figure 3-b). At instability threshold, we 
systematically recover two symmetrical points (/, q) and 
(/, —q) corresponding to the two helicoidal modes form- 
ing the braid pattern. Further away from the thresh- 
old, the 2D-Fourier transforms are more complex and 
indicate other peaks, but in the following we only fo- 
cus on the most energetic of them. As suggested by the 
analytical study of Shalybkov and Riidiger [13], we did 
not observe any systematic variation of the wavenumber 
k = 27r q {Ri{Ro — Ri)Y^'^: over all our experiments, we 
measure a mean value of 10.6 with a standard variation 
of 1.5 in close agreement with analytical and numerical 
results [3 • This value corresponds to a wavelength 
equal to 1.17 time the gap, smaller than the standard 
value in the absence of stratification (i.e. twice the gap) 
and confirms the reduction of the vertical extension of the 
structures because of stratification. Figure 4-a) presents 
the evolution of the experimentally determined pulsation 
uo = 27r/ / Qi divided by the mean normalized angular ve- 
locity, as a function of /i. As can be observed, the data 
points are clustered along horizontal lines corresponding 
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FIG. 4: (a) Evolution of the pulsation oo of the prevalent heli- 
coidal mode divided by the mean angular pulsation < Q(r) > 
for Fr = 0.5 obtained for N ranging from 0.85rad/s to 
S.lrad/s: all ratios align with integer values as expected from 
the azimuthal periodicity in the cylindrical geometry. These 
are highlighted by horizontal lines and correspond to the az- 
imuthal wavenumber m (from 1 to 4) of the selected mode. 
These values are used in (b) to rescale oo. All experimental 
data points correctly ali gns with the analytical predictions 
(dotted lines) of the SRI [11. 



to integer values. Indeed, extending the simple advection 
argument given by Caton et al. [6,] to /i > 0, one can esti- 
mate that the helicoidal patterns are simply advected at 
the mean angular velocity < fl{r) > of the fluid between 
the cylinders. Hence the ratio uj/< rt{r) > actually indi- 
cates the azimuthal wavenumber m of the selected mode, 
ranging from m = 1 to m = 4. We confirm these conclu- 
sions by directly measuring the azimuthal wavenumber 
at threshold when there is a single helix in the flow. In 
figure 4-b), we then compare our experimental results 
using the previously determined m with the numerical 
results of [I^. Again, perfect agreement is found. 
In conclusion, we have extended previous experimental 
studies of stratified Taylor- Couette flow by a systematic 
exploration of the stable range according to the Rayleigh 
criterion. We have observed the new Strato-rotational 
Instability in the form of non-axisymmetric oscillating 
modes. Our results validate the theoretical results of 
Molemaker et al. [H], Yavneh et al. 12], Shalybkov 
and Riidiger [l^ and Dubrulle et al. [14]. In particu- 
lar, instability thresholds, axial and azimuthal wave num- 
bers and temporal frequencies of the helicoidal travelling 
modes (m 7^ 0) have been measured for a fixed Froude 
number Fr = 0.5, an aspect ratio r] = 0.80 and for dif- 
ferent Reynolds numbers. Comparisons with theoretical 
predictions [l3] are excellent. Moreover, let us empha- 
sis that in the inviscid limit, the SRI threshold tends 
to the value /i = 7^, contrary to ji = rf given by the 
Rayleigh criterion valid in the unstratified case. As al- 
ready remarked by Shalybkov and Riidiger 0] , this im- 
plies that stratified Keplerian flows such as those found 
in accretion disks and characterized by 17 (r) ^ (i.e. 
l-^Kepier = rf^'^) are stable following the Rayleigh crite- 
rion (i.e. fiKepier > v"^)^ but cau be destabilized by the 

SRI as /J^Kepler < V' 
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